Scaling of the magneto conductivity of silicon MOSFET's: 
evidence for a quantum phase transition in two dimensions. 
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For a broad range of electron densities n and temperatures T, the in-plane magnetoconductivity 
of the two-dimensional system of electrons in silicon MOSFET's can be scaled onto a universal curve 
with a single parameter H„(n, T), where Ha obeys the empirical relation Htr = A{n)[A{n)'^ + T^Y'^ . 
The characteristic energy fcsA associated with the magnetic field dependence of the conductivity 
decreases with decreasing density, and extrapolates to at a critical density no, signaling the 
approach to a zero-temperature quantum phase transition. We show that H„ — AT for densities 
near no. 

PACS numbers: 71.30.-Hh, 73.40. Qv, 73.50.Jt 



Strongly interacting two-dimensional systems of elec- 
trons (or holes) have drawn intensive recent attention 
lyi due to their anomalous behavior: their resistance ex- 
hibits metallic temperature-dependence above a critical 
electron |g,|[ (or hole [Q-^) density Uc, raising the possi- 
bility of an unexpected metallic phase in two dimensions. 
An equally intriguing characteristic of these systems is 
their enormous response to magnetic fields applied in the 
plane of the electrons: the resistance increases dramati- 
cally with in-plane magnetic field and saturates to a new 
value above a characteristic magnetic field H^at on the 
order of several tesla 0-0. For high electron densities, 
measurements of Shubnikov-de Haas oscillations 11-1^] 
have established that the magnetic field Hgat is the field 
at which full polarization of the electrons is reached. A 
parallel magnetic field has been shown to suppress the 
metallic temperature-dependence fiH- Data obtained 
by Pudalov et al. fl^ and by Shashkin et al ||l3] indi- 
cate there is a substantial increase in the g-factor as the 
electron density is decreased toward Uc- These experi- 
mental findings all suggest that the behavior of the spins 
is key to understanding the enigmatic behavior of dilute, 
strongly interacting systems in two dimensions. 

In this paper we report measurements of the 
temperature-dependence and density-dependence of the 
in-plane magnetoconductivity of silicon metal-oxide- 
semiconductor field eflFect transistors (MOSFET's). For 
a broad range of electron densities and temperatures, we 
show that all data for the magnetoconductance can be 
collapsed onto a single universal curve using a single pa- 
rameter H^ which obeys an empirical relation given by 
H^{n, T) = A(n)[A(n)2 + T2]1/2, The characteristic en- 
ergy fcsA associated with the response to magnetic field 
is found to decrease with decreasing electron density, and 
to exhibit critical behavior, extrapolating to at a den- 
sity no near the critical density n^ for the zero-field metal- 



insulator transition. Hu = AT for densities near ng, 
so that the magnetoconductivity scales with H/T. Our 
results provide strong experimental evidence for a zero- 
temperature quantum phase transition at density uq. We 
suggest that this is a transition to a ferromagnetically or- 
dered state in two dimensions. 

Measurements were taken on three silicon MOSFETs: 
the mobility fi « 30,000 Y/icm^s) at 0.3 K for sam- 
ple #1 and « 20, 000 Y/{cm^s) for samples #2 and #3. 
Data were obtained on samples with split-gate geome- 
try at City College to 12 T and at the National Magnetic 
Field Laboratory in fields up to 20 T using standard four- 
terminal AC techniques described elsewhere [|2| . 

Fig. 1 (a) shows the longitudinal conductivity cr^x as 
a function of magnetic field i/|| applied parallel to the 
plane of a silicon MOSFET for different electron densi- 
ties 71 > Tic ~ 0.85 X 10^^ cm^^ (ric is the critical density 
for the zero-field metal-insulator transition). In agree- 
ment with earlier results, the conductivity decreases (re- 
sistivity increases) dramatically with increasing magnetic 
field and saturates to a value a{H = oo) that is almost 
independent of magnetic field in fields H > H„. 

Attempts 



10,15- 17 1 to obtain a collapse of the mag- 
netoresistance onto a single scaled curve have generally 
yielded scaling at either low or high magnetic field, but 
not over the entire field range. We now report a method 
that allows a full mapping of the magnetoconductiv- 
ity (rather than the magnetoresistivity) onto a single 
universal curve over the entire range of magnetic field 
for all data obtained at different densities n and tem- 
peratures r, as follows. We separate the conductivity 
into a field-dependent contribution [a — a{H = oo)], 
which was shown ||ll|,|l^ to reflect the degree of spin 
polarization of the electrons, and a contribution that 
is independent of magnetic field, a{H = oo). We find 
that the field-dependent contribution to the conductiv- 



ity, a{H = 0) - 
ajH = 0) - (j{H 
M, so that: 



<y{H), normalized to its full value, 
= cxd) is a universal function of H/ H^, 



cj{H = 0) -(t{H) 
a{H = 0) - (j{H = cx)) 



= F{H/H„) 



(1) 



using a single scaling parameter Ha-{n, T) which is a func- 
tion of electron density n and temperature T. Applied to 
the magnetoconductance curves shown in Fig. 1 (a) for 
different electron densities, the above scaling procedure 
yields the data collapse shown in Fig. 1 (b). 



H (tesla) 
10 




n 


(10" cm-') 






1.59 






1.88 






2.15 
2.43 
2.71 






2.99 




: 


3.34 
3.68 
4.38 



H/H 

a 

FIG. 1. (a) Conductivity versus in-plane magnetic field 
at different electron densities n in units 10^^ cm~^, as 
labeled. Data are shown for sample #3 at T = 100 
mK. (b) Data collapse obtained by applying the scal- 
ing procedure, Eq.(l), to the curves of Fig.l (a). 



Similar scaling holds for curves obtained at different 
temperatures. This is demonstrated in Fig. 2, which 
shows the scaled magnetoconductance for sample #2 
taken at a fixed density and different temperatures. 

There are small departures from a full data collapse 
with temperature and as the density is varied. Possi- 
ble reasons include: (a) the fact that the magnetocon- 
ductance is not strictly constant at high magnetic field, 
due possibly to orbital effects or disorder; (b) interfer- 
ence (weak localization) effects may become important at 



very high densities; (c) one expects scaling to hold only 
over a restricted range near a transition. Nevertheless, 
within an accuracy of about 5%, the magnetoconduc- 
tance curves for all electron densities and temperatures 
presented here can be mapped onto a universal function 
using a single parameter H/ H„. 




FIG. 2. Data collapse obtained by applying the scaling 
procedure, Eq.(l), to the in-plane magnetoconductivity 
of the 2D electrons at different temperatures for elec- 
tron density n^ = 0.94 x 10"'^^ cm~^. The inset shows 
the conductivity at different temperatures as a function 
of magnetic field. Data are shown for sample#2. 



Having thus demonstrated that a one-parameter scal- 
ing description of the magnetoconductance holds for tem- 
peratures up to 1.6 K over a broad range of electron den- 
sities up to 4 ric, we report below the results of a detailed 
experimental investigation of the behavior of the scaling 
parameter, iJo-, as a function of temperature and elec- 
tron density. In the analysis that follows, the values of 
He were determined by the intersection of the two dot- 
ted lines shown in Fig. 1 (b); this intersection signals 
a crossover from strong field-dependence at low fields, 
H < H„ , to a conductivity that is almost independent of 
magnetic field aX H > H„. 

Fig. 3 (a) shows the scaling parameter H^ plotted as 
a function of temperature for different electron densi- 
ties above the critical density for the zero-field metal- 
insulator transition. Tic. (Data are shown for sample^2, 
for which ric = 0.85 x 10^^ cm~^.) The scaling parameter 
becomes smaller as the electron density is reduced; for a 
given density, H„ decreases as the temperature decreases 
and approaches a value that is independent of 
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FIG. 3. (a) Ha as a function of temperature for differ- 
ent electron densities; the solid lines are fits to Eq.(2). 
(b) The parameter A versus electron density; the solid 
line is a fit to the expression A = Ao(7t. — no)". Data 
are shown for sample #2. 



temperature, Ha{T = 0). As the density is reduced to- 
ward Tie, the temperature-dependence of H^ becomes 
stronger and its low-temperature asymptotic value be- 
comes smaller. Note that for electron densities below 
1.36 X 10"'^^ cm~^, Ha is approximately linear with tem- 
perature at high T. The behavior of the scaling param- 
eter Hair) can be approximated by an empirical fitting 
function: 



10% of the critical density Uc ~ 0.85 x 10^^ cm~^ for 
the metal-insulator transition obtained from zero-field 
transport measurements. Additional measurements are 
required to determine whether no and ric are the same or 
different densities. 

The parameter A represents an energy fc^A and a cor- 
relation time th ~ ^/fcsA). For high densities and low 
temperature, T < A ^ I^/th, H^ is determined by A 
and the system is in the zero temperature limit. Near 
no the measuring temperature T > A '^ ^/th, the field 
Ha is dominated by thermal effects and is not in the 
T = limit. At n = no, the energy fc^A vanishes and 
the correlation time th diverges; the parameter Ha — > 
as T ^ 0; the system exhibits critical behavior, signal- 
ing the approach to a new phase in the limit T=0 at a 
critical density no M. 
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FIG. 4. The magnetoconductivity as a function of H/T 
of sample #1 at a density n = 0.82 x 10^^ cm^^, just 
below the critical concentration no The inset shows the 
conductivity as a function of magnetic field at different 
temperatures. 



i?,(n,T) = A(n)[A(n)2+T2]i/2 



(2) 



The solid lines in Fig. 3 (a) are fits to this expression 
using A{n) and A(n) as fitting parameters. As can be 
inferred from the slopes of the curves of Fig. 3(a), the 
parameter A{n) is constant over most of the range and 
then increases measurably at lower densities (less than 
20% within the measured range). As shown in Fig. 3 
(b), the parameter A decreases with decreasing density 
and extrapolates to zero at a value n = no. Similar 
behavior was found for sample#l (not shown). Within 
the experimental uncertainty of our measurements and 
of the analysis which yields A, a fit to the critical form, 
A = Ao(n — no)", yields a value for no that is within 



As discussed earlier, a number of experiments have 
shown that electron spins play an important role in the 
response of dilute 2D systems to an in-plane magnetic 
field. For high electron densities, Shubnikov-de Haas 
measurements |ll[-|l3| have shown directly that the un- 
usually strong field-dependence of the conductivity is as- 
sociated with the polarization of the electron spins by the 
magnetic field. The fact that the same scaling, Eq. (1), 
that applies at these high densities remains valid at low 
electron densities implies that spin polarization continues 
to be a major factor in the observed magnetoconductiv- 
ity as no is approached. The parameter Ha, the field that 
characterizes the system's response to in-plane magnetic 



field, vanishes at uq. Thus, the behavior of the magne- 
toconductivity as A —> in the vicinity of ng indicates 
critical behavior for the spin susceptibility. 

At no, A = 0, the scaling parameter H^ = AT, and the 
magnetoconductance scales with H/T. This is illustrated 
explicitly in Fig. 4, where the magnctoconductivity of 
sample #1 is shown as a function of H/T for a density 
near rig. Magnetization that scales with H/T could be 
associated with localized independent spins in a param- 
agnetic insulator, which would imply that the transition 
at no is to a localized phase. We note, however, that 
at densities near no the resistivity changes from weakly 
insulating to metallic temperature dependence, indicat- 
ing that the localization length is larger than the average 
distance between the spins. The spins can therefore not 
be considered as independent. 

It is possible that the critical behavior we report in this 
paper is due to a magnetic instability in the 2D Fermi sys- 
tem at density n = uq. Support for this scenario is pro- 
vided by the unusually strong temperature-dependence of 
Hct for metallic densities near ng and the increase found 
in the value of m*g* as the electron density is decreased 
||ll| , |l5|Jl6| , both of which are characteristic features of 
a metal near a ferromagnetic instability. As shown ex- 
perimentally by Okamoto et al. [ [ll| and Vitkalov et al. 
|l2| , |l3| , at high electron densities the ground state of the 
system corresponds to complete polarization of the elec- 
tron spins in high magnetic fields H > H„. The mag- 
netoconductance continues to scale, obeying Eq. (1), 
as the electron density is reduced toward no, implying 
there is no change in the character of the ground state at 
fields H > H„ as H^ decreases. This suggests that when 
H„ vanishes in the limit of zero temperature at n = no, 
parallel alignment survives at least on a scale of the lo- 
calization length Lioc- On this basis we suggest that a 
transition occurs at no due to a ferromagnetic instability. 
This could be associated with long range order or local 
ordering of ferromagnetic domains or droplets. 

In summary, we have shown that the in-plane magne- 
toconductivity of the dilute 2D electron system in silicon 
MOSFET's can be scaled with a single scaling parame- 
ter Her over a substantial range of electron densities and 
temperatures. The parameter H„{n,T) that character- 
izes the response of the system to a magnetic field obeys 
an empirical relation H^r = A(n)[A(n)^-|-T^]^/^. The en- 
ergy fc^A is found to be to zero at a critical density no, 
indicating the quantum critical regime of the 2D dilute 
electron system at n = no. These findings provide strong 
evidence for a zero-temperature quantum phase transi- 
tion at density no near (but not necessarily the same as) 
the density nc for the apparent metal-insulator transi- 
tion. 
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